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Abstract 

We report a numerical study of transport properties of a quantum dot with 

superconducting leads. We introduce a general phenomenological model of 

quantum dot transport, in which electron tunnel rates are computed within 

the Fermi's Golden Rule approach. The low temperature current-voltage (I- 

V) characteristics are in qualitative agreement with experimental observations 

of Ralph et al. [Phys. Rev. Lett., 74, 3241 (1995)]. At higher temperatures, 

our results reveal new effects due to the thermal excitation of quasiparticles 

in the leads as well as the thermal population of excited quantum levels in the 

dot. We also study the photon-assisted tunneling phenomena in our system 

and point out its potential for millimeter wave applications. 
PACS numbers: 73.40.Gk, 73.20Dx, 73.50Pz, 75.40.Mg 
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Due to the advancement of modern lithographic technologies, electrons can now be con- 
fined to a spatial region so small that not only the Coulomb charging energy Eq, but 
also the spacing of the quantum mechanical energy levels e, become accessible to experi- 
ments conducted at low temperatures. While the single electron charging effect has been 
observed in both metal tunnel junction systems^ and semiconductor quantum dots,i the 
discrete quantum levels have been studied mostly in semiconductor quantum dotsB! Re- 
cently, Ralph, Black and Tinkham conducted experiments in which they directly measured 
tunneling transport through discrete quantum levels in a small metallic particleJ! 

In this paper we present numerical simulation results based on the system studied by 
Ralph et a/.i Their samples were made of a tiny aluminum particle (with diameter < 10 
nm) connected to aluminum leads through tunnel junctions on both sides. A third electrode 
can also be added to act as a gate. Since the experiment was conducted at temperatures 
well below the superconducting transition temperature of aluminum (To = 1.21 K for their 
sample), the leads are nominally superconducting unless a magnetic field is applied. The 
dot itself can have discrete quantum levels with spacing either smaller or larger than the 
superconducting energy gap of aluminum. In this paper, we will assume that the quantum 
dot is a normal metal particle, and the leads can be either superconducting or normal. 

We introduce a general phenomenological model of transport through a quantum dot. In 
this model, we assume that the quantum dot is weakly coupled to the two leads by tunnel 
barriers. When an appropriate bias voltage V is applied to the leads, an electron can tunnel 
across one barrier into the dot and subsequently tunnel out through the second barrier. 
According to general tunneling theory,!'! the tunneling rate across a barrier from side "a" 
to side "b", can be evaluated using the Fermi's Golden Rule, 

r a ^b(/i a , = y f-oc \T ab \ 2 Af a (E - fl a )Af b (E - fl b ) 

xf(E-ii a )[l-f(E-to)]dE, (1) 

where T ab is the phenomenological tunneling matrix element, and f(x) = l/[l + exp(x//cBT)] 
is the Fermi function.00 ]\f a {E) and M b {E) are the density of states, and /i a and jj b are the 



chemical potentials, on their corresponding sides. For our system, to compute the tunneling 
rate from one of the leads to the dot, we take the BCS quasiparticle density of stateslll in 
the lead and assume that the dot itself has an evenly spaced (with spacing e) discrete level 
spectrum. Thus the density of states in the superconducting leads is, 

. if \E\ < A; 

Af s (E)=Af N { . ' ' (2) 



\E\ 



VE 2 -A 2 



if \E\ > A, 



and the density of states in the dot is, 

N D {E)=N N eY,KE-ne). (3) 

n 

Here, 2A is the superconducting energy gap of the leads, and Ajv is the density of states in 
the bulk normal metal that comprises both the leads and the dot. 

Once we have the tunneling rates, we can compute the tunneling current through the dot 
using a master equation approach.! At steady state, the current that flows into the dot from 
one lead should balance the current that flows out of the dot into the other lead. Therefore 
the steady state current is given by, 

I = [^L-*D{fiL, Vk)Pk-l - r D -> L (fJLk, fiL)Pk] , (4) 
k 

where //£ and are, respectively, the chemical potentials of the left lead and the quantum 
dot when the dot contains k electrons. The probability of the dot having k electrons is 
denoted by Pk, and it is determined from a steady state rate equation. iH In Eq. (f|), the 
first term in the summation is the contribution to the current from the k-ih electron on the 
dot which tunnels in from the left lead, given that the dot already has k — 1 electrons on it. 
The second term corresponds to the reverse process. The net total current is the sum over 
all possible k values, which depends on the applied voltage and the temperature. A more 
detailed account of our model and its application to various quantum-dot configurations are 
discussed elsewhere.il! We now discuss results concerning the particular system of a quantum 
dot with superconducting leads. 
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In Fig. [l|, we show a typical low temperature current- volt age (I-V) characteristic of our 
system. Here the temperature k%T = 0.02E C (E c = Ec + e is the spacing between chemical 
potential levels, and Eq = e 2 /Cs is the charging energy), the superconducting energy gap 
2A = 0.3E C and the quantum energy level spacing in the dot e = 0.2E C . When the leads 
are superconducting (solid curve), the I-V curve consists of a series of sharp peaks spaced 
e apart .0 This is in contrast with the I-V curve of the same dot with normal metal leads 
(dashed curve), which has only gentle steps with the same spacing e. 

The shape of these I-V curves and the various tunneling processes involved can be un- 
derstood using the illustration shown in the inset, which is a sketch of the energy spectra of 
the dot and the leads. In this sketch, the narrow vertical lines represent the two tunneling 
barriers. The quasiparticle density of states for the superconducting leads are shown; the 
shaded region indicates occupied states and the unshaded, unoccupied states. The chemical 
potentials of the superconducting leads are shown as dotted lines in the center of the su- 
perconducting gap. The energy spectrum for the dot includes both the excitation spectrum 
(dotted horizontal lines) and the addition spectrum (heavy horizontal lines). The excitation 
spectrum has spacing e and can be populated by, for example, thermal excitation of electrons 
within the dot without changing the electron number N. The addition spectrum, on the 
other hand, has spacing /ijv+i — /ijv = Eq, since it includes the Coulomb charging energy 
that one has to overcome in order to add an electron (N —>■ N + 1) to the dot by means 
of tunneling. Assume at zero temperature and zero voltage that there are N electrons on 
the dot and yP L and /i R are midway between /j, N and Hn+i- A dc voltage, eV — /i£ — fj, R is 
applied to the left lead while the right lead remains at fixed voltage. In order for an electron 
to tunnel into the dot from the left lead, one needs to raise //£,, such that \il — A > //jv+i, 
or V > V t = E c /2 + A. As soon as this threshold voltage, V t , is reached, an electron on the 
occupied states of the left lead can tunnel into the dot at level /xjv+i an d subsequently tunnel 
out to the right lead, resulting in a current flow which reflects the density of states of the 
left leadS As the voltage is increased further, the first excited state of the (N + l)-electron 
system (first dotted line above fiN+i) is reached. At this point, an electron can tunnel into 
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the dot either through level hn+i or the first excited state above it, before tunneling out into 
the right lead. This accounts for the increased current at a voltage e above the threshold 
voltage and having the shape again reflecting the density of states of the lead. We also see 
similar current rises corresponding to the second and higher excited states. However, as long 
as hl — A < Hn+2, the conduction cycle is always iV — > iV + 1 — > N . When the voltage is 
high enough, such that /i^ — A > hn+2 (eV > eV t + Eq), one extra electron can tunnel into 
the dot (through fiN+2 or the excited states above it) before the (N + l)-th electron can 
tunnel out. Now the electron number on the dot is a random choice among N, N + 1 and 
N + 2 (with corresponding probabilities P/v, P/v+i and Pn+2)- This corresponds to a larger 
increase in the current at the voltage V « 1.7. As the voltage is increased further, N will 
fluctuate among more and more possible configurations, and eventually the I-V will become 
essentially linear due to these fluctuations. 

The different current response with the superconducting and normal leads is simply 
the reflection of their respective density of states. The superconducting I-V has the sharp 
current rises and a higher threshold voltage because its density of state contains a gap and 
is singular at energies \E\ = A, both of which disappear in the normal leads where A = 0. 
Quantum dots with discrete levels and normal leads have been studied previously by Averin 
et alE 

The I-V curves become more complicated at higher temperatures due to thermal exci- 
tation of quasiparticles in the leads and electron population of the quantum dot excited 
states. Noting that there are many energy scales in our system and the fact that the super- 
conducting energy gap decreases with increasing temperature ,@ we will concentrate in the 
regime, Ec 3> k B T ~ [A(T),e]. This parameter range is clearly relevant to the experiment 
in Ref. |, where E c « 12 meV, e « 0.5 meV, A(0) = 0.18 and k B T c « 0.1 meV. 

In Fig. |2|(a) and (b), we show two I-V curves, both at relatively high temperature, such 
that in (a), 2A(T) > e, while in (b), 2A(T) < e. One of the intriguing features of these I-V 
curves is that in Fig. |2|(a) the narrow regions where the current is significantly suppressed 
do not correspond to the superconducting energy gap 2A(T). In fact, 2A(T) is much larger 
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as indicated in the figure. In Fig. 0(b), however, the relatively wide regions of current 
suppression do correspond to 2A(T). This seemingly contradictory phenomena has a simple 
explanation using the energy diagrams shown in the insets of Fig. 0, which are the finite 
temperature versions of the diagram that we used previously in discussing Fig. |I[ 

First of all, since [k B T/2A(T), ksT/e] ~ 1 (see the caption of Fig. 0), we expect that 
there will be many quasiparticles that are thermally excited across the gap. Also, quite a 
few quantum levels in the dot near the chemical potential levels will be partially occupied 
due to thermal excitations within the dot. Therefore many quantum levels, both below 
and above the chemical potential positions of the dot, will become available for tunneling. 
However, due to conservation of energy, a level can contribute to the current only if there 
are electrons in the left lead that have the same energy. Since there are no electron states 
within the gap, whenever a quantum level in the dot lies inside the gap of the left lead, 
that level will be excluded from conduction. When we sweep the voltage, if 2A(T) > e as 
in Fig. 0(a), there will always be either one or two quantum levels that fall inside the gap 
of the left lead. When the voltage is such that there are two levels lined up inside the gap 
and therefore excluded from conduction, the current will naturally drop in comparison with 
other voltage values where only one level is being excluded. This will explain Fig. 0(a). The 
situation in Fig. 0(b) is slightly different. Since we now have 2A(T) < e, there can only be 
one or no level being excluded from conduction. When a level falls inside the gap, it has 
to traverse the entire gap region before it can join the tunneling process again. Thus the 
current suppression regions in Fig. 0(b) have the same widths as the gap and are centered 
about every discrete level. 

Since the conduction in our system is carried out by tunneling processes, we expect 
there might be interesting ac properties to be explored when microwave radiation is coupled 
to the system. In particular, the photon-assisted tunneling phenomena, observed in large 
area superconductor-insulator-superconductor (SIS) tunnel junctions,0'El and semiconduc- 
tor quantum dotsJ^E should also manifest itself in some fashion in the present system. We 
treat our system in the presence of microwave radiation using the Tien-Gordon formalismjil 
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originally proposed to explain the photon-assisted tunneling phenomena observed in large 
area SIS junctions.0 Recently, the Tien-Gordon theory has also been applied to photon- 
assisted tunneling in semiconductor quantum structures.El'0 

We will treat the externally applied RF signal as an oscillating voltage, V r f cos(2irut), 
applied equally across each of the tunnel barriers. Then the Tien-Gordon theory leads to 
the following simple modification of the tunneling rates,0 

oo 

r(// -^ 6 ,a)= ^n( a ) r oOa- fi b + nhu) (5) 

n=— oo 

where a = eV r f j (hu), J n (x) are the Bessel functions of the first kind, and To is the tunneling 
rate in the absence of RF radiation, given in Eq. ([!]). 

In Fig. |3], we show I-V curves of our quantum dot with superconducting leads in the 
presence of RF radiation with frequency hu/e = 2/3 and two different power levels, along 
with the dark I-V (a = 0). Due to the periodic spikes in the dark I-V, its photo-response 
is considerably more complicated than that of the conventional SIS tunnel junctions, whose 
dark I-V has a single discontinuous jump at 2A.il For example, Fig. |3| illustrates that if the 
ratio hu/e = p/q is a rational number, then the I-V, in sufficiently strong RF signal, will 
have periodic spikes with period hu/p, not hu. Even more complicated behavior is observed 
when e and hu are incommensurate. We wish to point out that this effect, if confirmed by 
experiment, might find applications in millimeter wave detection and mixing schemes. In the 
weak signal limit, a < 1, one expects a single peak below the dark I-V threshold voltage V t , 
which will give information about the frequency as well as the amplitude of the RF signal. 
We found the sharpness of the current peaks in the superconducting-lead system makes it 
much more robust against thermal noise in comparison with its normal-lead counterpart. 

In summary, we have carried out numerical simulations of single electron transport 
through a quantum dot with superconducting leads. Our current- volt age characteristics 
calculated at low temperature is in good qualitative agreement with the experimental obser- 
vation of Ralph et a/.i In addition our analysis show that at higher temperatures, thermal 
excitation of quasiparticles in the leads and thermal population of the excited quantum levels 



within the quantum dot should lead to interesting changes in the I-V curves. We also predict 
that when RF radiation is coupled to the system, the photon-assisted tunneling phenomena 
should manifest itself by producing extra periodic structures in the I-V curves, which might 
be useful in the millimeter wave detector/mixer applications. Due to the presence of many 
different characteristic energy scales, the rich dynamical properties of this system demands 
more exploration. 
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FIGURES 

FIG. 1. Low temperature I-V characteristics of a quantum dot with superconducting leads 
(solid curve) and normal metal leads (dashed curve). The temperature is k B T = 0.02Eq, and the 
superconducting energy gap in the leads is 2A = 0.3Eq. The quantum level spacing is, e = 0.2Eq. 
The inset is a sketch of the energy spectra in the leads and the dot. Note the quantum dot energy 
spectrum includes the excitation spectrum (with spacing e), and addition spectrum (with spacing 
E* c = e 2 /C + e). 

FIG. 2. I-V characteristics at higher temperatures for (a) 2A(T) > e and (b) 2A(T) < e. The 
relevant energy ratios are: (a) T/T c = 0.815, k B T/E* c = 0.07, k B T/2A(T) = 0.311, k B T/e = 0.35, 
2A(T)/e = 1.14; (b) T/T c = 0.931, k B T/E* c = 0.08, k B T/2A(T) = 0.576, k B T/e = 0.40, 
2A(T)/e = 0.70. The vertical dotted lines are drawn at the positions of the quantum levels (e 
apart). 

FIG. 3. Photon assisted tunneling in quantum dots with superconducting leads. The dark I-V 
(solid curve) has the same parameter values as Fig. [l|(a). The dotted curve and the dashed curve, 
respectively correspond to power levels, a = eV r f /hf = 2.5, and a = 0.5, with hv/e = 2/3 for both 
curves. 
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Whan and Orlando Fig.2 
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